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Abstract. A search for pair-produced massive coloured scalar particles decaying to a four-jet final state is 
performed by the ATLAS experiment at the LHC in proton-proton collisions at yfs = 7 TeV. The analysed 
data sample corresponds to an integrated luminosity of 4.6 fb -1 . No deviation from the Standard Model is 
observed in the invariant mass spectrum of the two-jet pairs. A limit on the scalar gluon pair production 
cross-section of 70 pb (10 pb) is obtained at the 95% confidence level for a scalar gluon mass of 150 GeV 
(350 GeV). Interpreting these results as mass limits on scalar gluons, masses ranging from 150 GeV to 
287 GeV are excluded at the 95% confidence level. 

PACS. 12.60.Jv Supersymmetric models, - 14.80.Ly Supersymmetric partners of known particles 



Massive coloured scalar particles that decay into glu- 
ons are predicted in several extensions of the Standard 
Model (SM). The most prominent examples are the scalar 
partners of a Dirac gluino called scalar gluons ('sgluons') 
in extended supersymmetric models like the Af=l / Af=2 
hybrid model [Tj2] or the R-symmetric MSSM [3p]. These 
particles are also present in compositeness models [5j|9] 
where they are called hyperpions. While single produc- 
tion of sgluons is possible, the production cross-section 
depends strongly on the masses of the supersymmetric 
particles and, in typical supersymmetric scenarios, is of 
the same order as the pair production cross-section. On 
the contrary the pair production cross-section does not, 
at leading order, depend on supersymmetric parameters 
except for the sgluon mass. Since the sgluon has posi- 
tive R-parity [10| and since the sgluon coupling to quark- 
antiquark pairs is suppressed by the quark mass, light 
sgluons, i.e. sgluons with masses of the order of 100 GeV, 
are expected to decay to two gluons with a branching ra- 
tio close to one [2j[4] . Pair production of sgluons each de- 
caying to two gluons, leading to a four-jet final state, is 
therefore used as a benchmark process. ATLAS has pre- 
viously searched for signatures of these particles in the 
_1 recorded in 2010 11 , excluding at the 



dataset of 34 pb" 
95% confidence level (CL) sgluons with masses of 100 GeV 
to 180 GeV, with the exception of a mass window of 5 GeV 
around 140 GeV. The search described in this paper, us- 
ing data recorded in 2011, explores the mass region from 
150 GeV up to about 300 GeV. 

The strategy of the analysis is to first reconstruct the 
masses of the two sgluon candidates. Then the distribu- 
tion of the average mass of the two sgluon candidates is 
analysed for evidence of a signal with a fit to the back- 
ground plus a signal template of variable strength. The 



background for this search is the SM multijet background, 
due to its large cross section. 

ATLAS is a multipurpose detector [i"2p3] with nearly 
4-7T coverage in solid angle. The inner aetector, consist- 
ing of silicon pixel and microstrip detectors as well as a 
transition radiation tracker, is immersed in a 2 T axial 
magnetic field. In the pseudorapiditjr] region |ry| < 3.2, 
high-granularity lead/liquid- argon (LAr) electromagnetic 
(EM) sampling calorimeters are used. An iron/scintillator- 
tile calorimeter provides hadronic coverage for \rj\ < 1.7. 
The end-cap and forward regions, spanning 1.5 < \rj\ < 
4.9, are instrumented with LAr calorimeters for both EM 
and hadronic measurements. The calorimeters are sur- 
rounded by a muon spectrometer which consists of three 
large superconducting toroids, a system of precision track- 
ing chambers, and fast detectors for triggering. ATLAS 
uses a three-level trigger system. The first-level trigger is 
implemented in custom hardware, the other two trigger 
levels are implemented in software running on commer- 
cially available PC farms. 

SM multijet Monte Carlo (MC) samples are used to 
develop a data-driven method for the background es- 
timation. To incorporate detector effects, these events 
are passed through a detailed simulation of the ATLAS 
detector 



[14] 



based on Geant4 15 . Alpgen 16 SM 



multijet samples are generated with the MLM matching 



1 ATLAS uses a right-handed coordinate system with its ori- 
gin at the nominal interaction point (IP) in the centre of the 
detector and the z-axis along the beam pipe. The x-axis points 
from the IP to the centre of the LHC ring, and the y-axis 
points upward. Cylindrical coordinates (r, 4>) are used in the 
transverse (x-y) plane, <f) being the azimuthal angle around 
the beam axis. The pseudorapidity is defined in terms of the 
polar angle 9 as r\ = — lntan(#/2). 
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scheme [17] and interfaced to Herwig [18] for the parton 



shower and fragmentation processes and to Jimmy 19 
for the simulation of the underlying event. The Alpgen 
samples are generated with the CTEQ6L1 parton dis- 
tribution functions (PDFs) [20] (underlying event tune 
AUET2-CTEQ6L1 [2l]). Pythia [22] SM multijet sam- 
ples are generated with the LO* MRST PDFs [23] (un- 
derlying event tune AUET2B LO** [2l]), The sgluon pair 
production differential cross-section T2J is implemented as 
an external process interfaced to Pythia. The decay of 
the sgluons as well as fragmentation and hadronisation are 
performed by Pythia. For the normalisation of the signal 



a next-to-leading-order (NLO) cross-section is used 24 



Signal samples of 40k to 120k events each were gener- 
ated with masses of 150, 200, 250, 300 and 350 GeV and 
passed through the fast simulation of the ATLAS detec- 
tor. Additionally the fast simulation was verified by pass- 
ing some samples through the detailed simulation which 
leads to results with good agreement. 
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Fig. 1. The reconstructed average mass distribution after all 
analysis cuts for the signal samples with m sg i uon =150, 250 and 
350 GeV. The curves are normalized to unit area. 



Jets are reconstructed using the anti-fcj jet cluster- 
with a radius parameter of 0.4. The 



25 



ing algorithm 

inpu ts t o the jet algorithm are three-dimensional clus- 
formed from energy deposits in the calorimeters. 
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ters 

The jets are calibrated using transverse momentum (px) 
and ry-dependent correction factors based on MC simula- 
tions and validated by test beam and collision data stud- 
Quality criteria are applied to reject jets produced 

Such jets are typically 



ICS 
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by non-collision backgrounds 
produced by hardware problems in the calorimeters, LHC 
beam-gas interactions or cosmic-ray induced showers. 

The analysis uses collision data collected in the year 
2011 at a centre-of-mass energy of y/s = 7TeV and cor- 



responding to an integrated luminosity of 4.6 fb . The 
data were recorded with a multijet trigger requiring at 
least four jets. Reflecting the trigger px thresholds, the 
trigger efficiency does not depend strongly on the px of 
the four highest-px jets for px > 80 GeV. To ensure full 
trigger efficiency for the analysis, the four highest-px jets 
in a selected event are required to be separated from each 
other by AR > 0.6 (where AR = y/{A4>) 2 + (Ar/) 2 with 
A<j) and Ar\ being the difference in azimuth and pseudora- 
pidity of the two jets). The resulting trigger efficiency of 
at least 99% obtained with simulated events was verified 
in data with the use of a single-electron trigger. 

A selected event must contain at least one recon- 
structed primary vertex with five or more associated 
tracks each having px > 400 MeV. At least four jets 
with pt > 80 GeV and \rj\ < 1.4 are required, since 
the signal is produced centrally, in contrast to the SM 
multijet background. These selection criteria together 
with the trigger requirements are referred to as prese- 
lection in the following. To improve the sensitivity of 
the analysis, the jet px threshold is defined as a func- 
tion of the probed sgluon mass. The of the fourth 
highest-px jet is required to be greather than the max- 
imum between 80 GeV and 0.3 times the sgluon mass plus 
30 GeV (p T (4*' l jet) > max(0.3 x m sgluon + 30, 80) GeV). 
The stringent px requirements select sgluons produced 
with a boost and hence lead to an accumulation of jet pairs 
from sgluon decays with AR w 1. Taking advantage of this 
property and to reconstruct the two sgluon candidates, 
the four highest-px jets in a preselected event are paired 
by minimising the quantity |Z\i? pa i r i — 1| + |ARpair2 — 1|. 
The event is rejected if, for the chosen combination, a jet 
pair has AR > 1.6. The reconstructed masses of the sglu- 
ons are denoted mi and 777,2. The reconstructed average 
mass is (mi + m 2 )/2. The scattering angle is defined as 
the angle between the direction of motion of the recon- 
structed sgluons in the rest frame of the four highest-px 
jets and the boost direction between the lab frame and 
the rest frame of the four highest-px jets. The magnitude 
of the cosine of the scattering angle | cos(6>* )| is required 
to be less than 0.5. In fact the SM multijet background 
is peaked in the forward region, reflecting ^-channel gluon 
exchange, while the signal is produced centrally due to 
the (mainly) s-channel production and scalar nature of 
the sgluon. Finally, to further improve the rejection of the 
SM multijet background, the relative difference between 
the two reconstructed masses \mi — m^/irni + 777.2) is re- 
quired to be less than 0.15. The selection efficiency of the 
signal is 0.6% for all simulated samples of the different 
sgluon masses. 

The reconstructed average mass distribution after all 
cuts is shown for sgluon signals with varying masses in 
Fig. [T] The natural width of the sgluon in this mass range 
is negligible and the width of the observable is entirely 
dominated by the instrumental mass resolution. As the 
mass of the sgluon increases, the requirement on the trans- 
verse momentum of the jets becomes relatively less strin- 
gent. An increase of 100 GeV in the sgluon mass leads to 
an increase of only 30 GeV in the requirement on the jet 
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transverse momenta. As a consequence the radiative tails 
to lower masses are more evident for higher masses in the 
figure, as they are less sculpted by the pr cut. 

After applying preselection and jet pairing, the pri- 
mary variables used in the analysis are compared to the 
Alpgen and Pythia multijct simulation. Backgrounds 
other than SM multijet events are estimated to be smaller 
than 1% of the total background sample and are thus ne- 
glected. Since the analysis requires at least four jets, Alp- 
gen is expected to give a better description than Pythia, 
which generates the third and fourth jets via a parton 
shower. 

The MC samples are normalised to the data after the 
preselection described above. The normalisation factor of 
1.25 obtained for Alpgen is compatible with the one 
(1.26) obtained in Ref. 11 . For Pythia the normalisa- 
tion factor is 0.75, close to the value of 0.65 [29] obtained 
with the 2010 data using a different tune (AMBT1). After 
this normalisation, for the transverse momentum of the 
4 th jet (Fig. p2la)), the separation of the two jets of the 
sgluon candidate with the highest transverse momentum 
(Fig.[2]jb)), the relative mass difference (Fig. 2jc)) and the 
cosine of the scattering angle (Fig. [2ld)), the agreement 
between the data and the MC simulations is, in general, 
at the 20% level. To reduce the dependence on simulation, 
the background is estimated from data, taking advantage 
of the kinematic properties of the sgluon signal. 

The main discriminating variable for the analysis is 
the reconstructed average mass. To determine the back- 
ground normalisation as well as the shape of the back- 
ground in the signal region, an ABCD method is used. 
The data sample is divided into one signal region (A) 
and three background-dominated regions (B, C and D). 
The variables used to define the four regions are | cos(6>*)| 
and mi — rri2\/(mi + 7712). The regions defined in Ta- 
ble [T are chosen as a compromise between the statistical 
significance of the signal in region A and the statistical 
uncertainty in the regions B, C and D which feeds into 
the uncertainty on the background prediction. 



Table 1. Definition of the four regions for the background 
determination. Region A is the signal region. 



Region 


|cos(0*)| 


mi — m 2 /(mi + m 2 ) 


A 


< 0.5 


< 0.15 


B 


> 0.5 


< 0.15 


C 


< 0.5 


> 0.15 


D 


> 0.5 


> 0.15 



The correlation between the two variables is less than 
0.1% in the four regions in the data and less than 1% in the 
Pythia samples, so the normalisation of the background 
in the region A is derived from the ratio of events in the 
control samples using N' A = Nb • Nq/N^. A closure test is 
performed with the Pythia and Alpgen MC samples and 
shows that N' A reproduces the actual number of events in 
region A, Na, within 2%. The difference is assigned as a 
systematic uncertainty on the background prediction. 



Table [2] summarises the results obtained in data, for 
the five sgluons mass hypotheses and corresponding signal 
regions, A, together with the corresponding background 
predictions. Good agreement is observed. The assumption 
of the ABCD method that the shapes of the reconstructed 
average mass distributions in regions A and B are the same 
was verified on the Monte Carlo samples. The last column 
gives the p-value obtained from a Kolmogorov-Smirnov 
test between the shapes of the reconstructed average mass 
distributions for data in regions A and B. Satisfactory 
p-values are found in this test on the data, which considers 
statistical uncertainties only. 

The result of the background estimation is shown in 
Fig. [3] for sgluon masses of 150, 250, 300 and 350 GeV. 
The data in region A are compared to the data in the 
control region B normalised using the ABCD method. The 
expected sgluon signal in region A is also shown. The ra- 
tio of background to data as well as the significance, in 
standard deviations, of the difference between the data 
and the prediction are shown as a function of the recon- 
structed average mass. The significance takes into account 
only statistical uncertainties. 

The average mass distribution of the data in region A 
is compared to the background-only prediction and to the 
background-plus-signal prediction using a binned likeli- 
hood test, which incorporates the signal contamination in 
the control regions and systematic uncertainties via nui- 
sance parameters. Systematic uncertainties affecting the 
simulated sgluon signal shapes are incorporated in the fit 
by varying signal templates. 

The systematic uncertainties on the acceptance, and 
the correlation assumed for each uncertainty source be- 
tween the four regions, are listed in Table [3] for a sgluon 
mass of 300 GeV. The uncertainty on the jet energy scale 
as well as the uncertainty on the jet energy resolution im- 
pacts both the signal shape and the acceptance. These un- 
certainties were measured with the complete 2010 dataset 
[30] to which an extra uncertainty for the higher pile-up in 
tEe 2011 run was added. The variation of the Pythia pa- 
rameters controlling initial state and final state radiation 
in a range consistent with experimental data 31 produces 



only a small effect on the reconstructed average mass dis- 
tribution. Therefore, this systematic uncertainty is taken 
into account for the signal acceptance. An uncertainty of 
1% is assigned to the trigger efficiency for sgluon signals. 
The uncertainty on the signal acceptance due to the signal 
MC statistical uncertainty is uncorrelated among the four 
regions. The acceptance uncertainty due to the PDFs is 
estimated using the independent CT10 20 error sets. The 

De 3.9% [32l[33l. 



uncertainty on the luminosity is taken to 

To probe for the presence of a signal, a fit with a 
freely varying signal strength parameter is performed for 
each mass hypothesis. No significant deviation from zero 
is found, and limits on sgluon production are derived. 
The profile likelihood ratio q„ |34| is used as test 



statistic, and exclusions are determined using the CL S ap 



proach 35 . Exclusion limits are computed using samples 
of pseudo-experiments generated taking into account all 
uncertainties and also contamination of the control re- 
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Fig. 2. The kinematic variables of the analysis are shown after applying the preselection and pairing the four highest-px jets: 
(a) is the transverse momentum of the fourth highest-px jet. (b) is the AR between the two jets of the reconstructed sgluon 
candidate with the highest transverse momentum jet. (c) is the relative mass difference, (d) is the cosine of the scattering angle 
in the four-jet centre-of-mass frame. The black histogram is the signal for a sgluon mass of 150 GeV normalised to the NLO 
cross-section. Data (dots) are compared to the Alpgen (triangles) and Pythia (rectangles) SM multijet samples where the MC 
samples are normalised to the data. The ratio data/MC is also shown separately for Alpgen and Pythia with its statistical 
uncertainty. 



gions by signal events. The normalisation and shape of 
the event distribution in region B are used, whereas for 
regions C and D only the normalisation is used. MC tem- 
plates are used to generate the shape of the signal in each 
pseudo-experiment. In each pseudo-experiment the statis- 
tical and systematic uncertainties are randomised, using 
Poisson and Gaussian distributions. 



Figure [4] shows the expected and observed 95% CL up- 
per bounds on the product of the scalar pair production 
cross-section and the branching ratio to gluons as a func- 
tion of the scalar mass. For a mass of 150 GeV (350 GeV), 
a limit of 70 pb (10 pb) on the scalar gluon pair produc- 
tion cross-section at the 95% CL is obtained. The solid 
line corresponds to the prediction of the sgluon pair pro- 
duction cross-section at NLO 24 , which is larger than 
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Fig. 3. The comparison of the data in the signal region with the background prediction is shown for: (a) m ag i uon = 150 GeV, (b) 
Ti a gi uon = 250 GeV, (c) m S gi uon = 300 GeV and (d) m sg i uon = 350 GeV. The points are the data in the signal region (region A). 
The plain histogram (red) is the expected signal in region A normalised to the NLO cross-section. The prediction of background 
in region A based upon the data in region B normalised using the ABCD method is shown as the rectangles which include 
the statistical uncertainty. The data/background ratio and the statistical significance of its difference from one, in standard 
deviations, are also shown in the lower panels. 



the leading order cross-section by a factor of about 1.6. 
The hatched band indicates the systematic uncertainty 
due to the choices of renormalisation and factorisation 
scales. Due to this recent NLO calculation, the previously 
unexcluded mass region around 140 GeV 



11 



is now ex- 



cluded by reinterpreting the limits obtained with the data 
recorded in 2010. For the analysis of the data recorded 
in 2011, sgluons with a mass from 150 GeV to 287 GeV 
are excluded. The endpoint of the mass limit is defined as 



the intersection of the cross-section limit with the NLO 
cross-section minus one standard deviation of the theory 
uncertainty. 

The dashed line is the prediction for the hyperpion 
cross-section of a compositeness model, obtained by rescal- 
ing the sgluon cross-section according to the ratios from 
Ref. 0. Since the ratios were calculated at leading or- 
der, this line should only be considered as an approximate 
indication of the excluded mass region. 
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Table 2. Comparison of the data in the signal region with the background prediction. The first column is the sgluon mass 
hypothesis, the second column is the corresponding minimum pt requirement on the four jets, followed by the number of 
observed data events in the signal region A (third column) and the number of predicted events using the ABCD method (fourth 
column), where the first uncertainty given is statistical and the second systematic. The last column gives the p-value obtained 
from a Kolmogorov-Smirnov test between the shapes of the reconstructed average mass distributions in regions A and B. Only 
statistical uncertainties are considered in this test. 

Sgluon Mass [GeV] pff m [GeV] Data ABCD prediction Shape p-value(A,B) 



150 80 102162 101100 ±800 ±2000 0.22 

200 90 55194 54500 ± 600 ± 1100 0.10 

250 105 23404 22500 ± 340 ± 500 0.28 

300 120 11082 10640 ± 230 ± 210 0.24 

350 135 5571 5330 ±180 ±110 0.70 



Table 3. The systematic uncertainties on the signal due to the jet energy scale (JES), jet energy resolution (JER), initial and 
final state radiation (ISR/FSR), the trigger efficiency (Trigger), the Monte Carlo signal statistics (MC Statistics), the choice 
of parton distribution functions (PDFs) and the integrated luminosity (Luminosity). The relative uncertainty of the signal 
acceptance is given for the four regions and for a sgluon mass of 300 GeV. The JES uncertainty is treated as asymmetric, 
corresponding to upward and downward fluctuations of the JES. For the JER uncertainty, only an upward fluctuation of the 
JER is considered, i.e. only a degradation of the JER, with respect to the nominal MC JER. The last column shows the expected 
correlation among the four regions. 



Systematic 


A 


B 


C 


D 


Correlation ABCD 


JES 


±10%,-10% 


±11%,-11% 


±11%,-13% 


±15%,-10% 


100% 


JER 


±0%,-2% 


±0%,-7% 


±0%,-l% 


±0%,-2% 


100% 


ISR/FSR 


±3.5%,-3.5% 


±3.5%,-3.5% 


+3.5%,-3.5% 


+3.596,-3.5% 


100% 


Trigger 


±1% 


±1% 


±1% 


±1% 


100% 


MC Statistics 


±4% 


±11% 


±5% 


±8% 


0% 


PDFs 


±3% 


±3% 


±3% 


±2% 


0% 


Luminosity 


±3.9% 


±3.9% 


±3.9% 


±3.9% 


100% 



with the data-driven background estimation. No evidence 
for new phenomena is found. Cross-section times branch- 
ing ratio limits as a function of the mass of the scalar par- 
ticle are derived. Interpreting the cross-section limit with 
the NLO cross-section calculation, sgluons with masses 
from 150 GeV to 287 GeV are excluded at 95% CL us- 
ing the data recorded in 2011 while the range 150 GeV to 
306 GeV is expected to be excluded in the absence of a sig- 
nal. Combining the exclusions obtained with the analyses 
of 2010 and 2011 data, sgluons with masses from 100 GeV 
to 287 GeV are excluded at the 95% CL. 
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